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Abstract 

Hybrid inverse problems are mathematical descriptions of coupled-physics (also 
called multi-waves) imaging modalities that aim to combine high resolution with 
high contrast. The solution of a high-resolution inverse problem, a first step that 
is not considered in this paper, provides internal information combining unknown 
parameters and solutions of differential equations. In several settings, the internal 
information and the differential equations may be described as a redundant system 
of nonlinear partial differential equations. 

We propose a framework to analyze the uniqueness and stability properties of 
such systems. We consider the case when the linearization of the redundant system 
is elliptic and with boundary conditions satisfying the Lopatinskii conditions. 
General theories of elliptic systems then allow us to construct a parametrix for 
such systems and derive optimal stability estimates. 

The injectivity of the nonlinear problem or its linearization is not guaranteed 
by the ellipticity condition. We revisit unique continuation principles, such as 
the Holmgren theorem and the uniqueness theorem of Calderon, in the context of 
redundant elliptic systems of equations. 

The theory is applied to the case of power density measurements, which are 
internal functionals of the form 7|Viip where 7 is an unknown parameter and u 
is the solution to the elliptic equation V • 7Vn = on a bounded domain with 
appropriate boundary conditions. 



1 Introduction 

A recent class of (mostly medical) imaging modalities, called hybrid, coupled-physics, 
or multi-wave modalities offers the possibility to reconstruct high-contrast parameters 
of interest with high resolution. High contrast is important to discriminate between, 
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say, healthy and non-healthy tissue. Resolution is important to detect anomalies at an 
early stage. 

Such hybrid modalities typically involve two steps. In the first step, not considered in 
this paper, a high resolution modality takes as an input measurements performed at the 
boundary of a domain of interest and provides as an output internal functionals of the 
parameters of interest and of specific solutions of underlying partial differential equations 
describing the probing (medical imaging) mechanism. This paper is concerned with the 
second step, involving quantitative reconstructions of parameters from knowledge of said 
internal functionals. For recent books and reviews on hybrid inverse problems, we refer 
the reader to, e.g., H il El EHl SSI |19] . 

Most practically used hybrid inverse problems involve internal functionals that are 
polynomials in the parameters of interest and the specific solutions mentioned above. 
Combined with the equations describing the latter solutions, we observe that all avail- 
able information represents a coupled, often redundant, system of nonlinear partial 
differential equations. 

In some instances, local algebraic manipulations allow us to solve such a system ex- 
plicitly. In the framework of functionals of solutions to second-order equations (and not 
of their derivatives), we refer the reader to, e.g., [HI [121 Ull CSl HHj. Such theories find 
applications in the quantitative step of the imaging modalities Photo-acoustic tomog- 
raphy, Thermo-acoustic tomography. Transient Elastography, and Magnetic Resonance 
Elastography; see also [5T1 [221 ISH [23 [ID [IS] • In the framework of functionals of the gra- 
dients of solutions, which find applications in Ultrasound Modulated tomography and 
in Current Density Imaging, we refer the reader to, e.g., O [3 [HI [131 IIHl IMl [301 ESI [37] . 

In many cases, explicit algebraic inversions may not be known or may not be appli- 
cable because not enough information is available. This paper proposes a framework to 
address several such problems when the linearization of the coupled system is elliptic. 
Hybrid inverse problems need not be elliptic; see the example of the 0-Laplacian in 
[71 [19] (also recalled below in section 12.6. ip or the Photo-acoustic problem as treated in, 
e-g-! [m HSj. However, when the number of internal functionals increases, the result- 
ing hybrid system becomes more redundant and hence more likely to be elliptic. We 
consider such a setting in section [21 We recall that elliptic systems augmented with 
boundary conditions that satisfy the Lopatinskii conditions admit left-parametrices. 
This follows from the theory of Agmon-Douglis-Nirenberg [21 [3] and the extensions to 
redundant systems by Solonnikov [34] . The existence of parametrices allows us to solve 
the linear problem up to possibly a finite dimensional space. Along with the construc- 
tion of a parametrix, elliptic regularity theory provides optimal stability results for the 
linearization of the hybrid inverse problem. A similar setting, which does not account 
for boundary conditions, is considered in ^T\ . 

The possible existence of a finite dimensional kernel for the linearized hybrid in- 
verse problem prevents us from determining whether the available internal functionals 
uniquely determine the coefficients of interest. Moreover, the dimension of the finite 
dimensional kernel is not stable with respect to small perturbations, which prevents us 
from analyzing the uniqueness and stability properties of the nonlinear hybrid prob- 
lem. A powerful methodology to obtain uniqueness results in the framework of elliptic 
systems of equations is the notion of unique continuation. In section [3l we revisit two 
classical notions of unique continuation. One is based on the Holmgren theorem, which 
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we generalize to the setting of redundant systems considered in this paper. The second 
one is based on the use of Carleman estimates as they are formulated in the Calderon's 
uniqueness theorem. See [TTJ [HI |26l [271 [2H1 [331 [iQl [SO] for references on these unique 
continuation results. Several extensions of these results, following the presentation in 
[10], are given in the setting of redundant systems in section [HI with proofs postponed 
to the appendix. 

Once a reasonable uniqueness result has been obtained for the linearization of the 
nonlinear hybrid inverse problem, several statements about uniqueness and iterative 
reconstruction procedures can be formulated for the nonlinear hybrid inverse problem. 
A constructive fixed point iteration method and a non-constructive local uniqueness 
result for the nonlinear problem are presented in section [H 

As an application of the conditions of ellipticity including boundary conditions and 
the conditions for unique continuation, we consider the case of power density internal 
functionals Hj{x) = ■y{x)\Vuj\'^{x), where Uj is the solution to V ■j'Vuj = on an open 
domain X C with boundary conditions Uj = fj on dX for 1 < j < J. For such a 
problem, we characterize the conditions under which the redundant problem is elliptic 
(for J = 2 in dimension n = 2 and J = 3 in higher dimension) and analyze cases in 
which a unique continuation principle (UCP) applies. 



2 Inverse Problems with local internal functionals. 

2.1 Systems of nonlinear partial differential equations 

Let 7 be a set of constitutive parameters in (linear or nonlinear, scalar or systems of) 
partial differential equations of the form 

£(7, Uj) =0 in X, Buj = fj on dX, (1) 

where £ is a polynomial in the derivatives of the solution uj and those of 7 on the open 
domain X C M", with uj augmented with boundary conditions on dX for 1 < j < J. 
Let us now assume knowledge of the functionals 

A^(7, Uj) = Hj in X, 1 < J < J. (2) 

where is a polynomial in the derivatives of the solution uj and those of 7. 

Several hybrid inverse problems may be recast in this general framework. More 
generally, we could have knowledge of functionals of the form A^(7,Mj,Mj) = Hij, or 
functionals Ai depending on more than two solutions Uj. We restrict ourselves to (|2|) 
to simplify notation. 

The above problem may thus be recast as a system of nonlinear partial differential 
equations for (7, {uj}): 

C{'y,Uj) = in X, Buj = fj on dX, 'i- < j < J 
M{^,Uj) = Hj inX, 1<3<J. 

The first relevant question for the inverse problem is whether the above system 
admits a unique solution. Since the solutions Uj are uniquely determined by knowledge 



3 



of 7, we are primarily interested in finding a unique solution to the parameters 7. 
However, looking at the coupled system of equations for (7, {uj}) is more revealing. 

The second question pertains to the stability properties of the reconstruction. Pro- 
vided that the solution to the inverse problem is unique, we wish to understand how 
perturbations in the information {Hj} propagates to the reconstruction of 7. 

The uniqueness and stability properties of the system depend on the number of 
acquired internal functionals J and on the way the medium was probed via the boundary 
conditions {fj}- Understanding how the uniqueness and stability properties are affected 
by changes in J and the boundary conditions {/j} is the third question we wish to 
answer. 

2.2 Linearization 

Some problems of the form (|3]) can directly be solved as non-linear systems. For instance, 
when £ is a linear second-order equation in uj and A^(7,Mj) = Uj the solution itself, 
the full non-linear problem is analyzed in [HI ITS] . 

For many problems in which direct reconstruction procedures may not readily be 
available, it is fruitful to analyze the linearization of (|3]). Neglecting boundary conditions 
at first, this yields 

dyC{'y,Uj)S'y + duC{'y,Uj)Suj = in X, l<j<^ 
dyM{^,Uj)6'-f + duM{'y,Uj)6uj = 6Hj in X, 1 < j < 

Note that the above differential operators may all be of different orders. With v = 
(57, {Suj}y, we may recast the above system as 

Av = 5, (5) 

for some appropriately defined source S. Let us assume that each uj is a scalar solution 
and that each Hj is also a scalar information. Then ^ is a system of differential operators 
of size 2 J X ( J + Af), where M is the number of scalar functions describing 7. Note that 
for J < M, the above system is under-determined. We consider here the case J > m. 
In several practical problems, J = M gives rise to a determined system A that is either 
not invertible, or invertible with poor stability properties. It is therefore also fruitful to 
consider the setting with J > M. 

2.3 Ellipticity 

In applications, C is often a linear, elliptic, operator. Adding the constraints (|2]), how- 
ever, may render the coupled system ([3]) or its linear version non-elliptic. One fruitful 
strategy to solve on the whole domain X (with appropriate boundary conditions) is 
therefore to "ellipticize" A, i.e., to find a number of constraints J sufficiently large so 
that A is elliptic, provided that such a J exists. 

What we mean by elliptic is defined as follows. For each x E X, A{x, D)ij is a 
polynomial in D = [dx^, ■ ■ ■ ,dx„) for 1 < i < 2J and 1 < j < J + M. We define 
the principal part Aq of A in the sense of Doughs and Nirenberg [23]. For each row 
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1 < z < 2 J of the system, we associate an integer Sj and for each column 1 < j < J + M 
of the system an integer tj. We normahze these integers by assuming that max(sj) = 0. 

We assume that Aij{x,D) is a polynomial in D of degree not greater than Si + tj. 
Then Ao^ij{x,D) is the part of the polynomial in Aij{x,D) of degree exactly equal to 

■Sj ~l~ tj. 

When all differential operators in (jlj) are of the same order t, then we may choose 
Si = and tj = t, in which case Ao{x, D) is composed of entries that are homogeneous 
polynomials of degree t in D. Many practical problems arise in forms in which the 
differential operators in (jlj) have different orders. 

We say that A is elliptic when the matrix Ao{x,^) = {Ao,ij{x,C,)} is full-rank (i.e., 
of rank J + M here) for all ^ G §"~^ the unit sphere and all x G X. 

Being full-rank is more likely when J is large, i.e., when A is over-determined. It 
is then useful to acquire redundant information Hj until ([5]) above is elliptic because 
elliptic systems enjoy more favorable stability estimates than non-elliptic systems. 

2.4 Lopatinskii boundary conditions 

Let us assume that we have been able to prove that A was a redundant elliptic system of 
equations. Then the system can be solved, up to possibly a finite dimensional subspace, 
when it is augmented by boundary conditions that satisfy the Lopatinskii criterion. This 
is defined as follows; see |44j . 
We consider the system 

Av = S inX, Bv = (l) ondX, (6) 

where B{x,D) is a Q x ( J + M) matrix with entries Bij{x,D) that are polynomial in 
D for each x G dX. We denote by % the order of By and by cij = maxj{bij — tj). 
Then Bo{x,D) is the principal part of B and consists of entries Bo^ij{x,D) defined as 
the polynomials of Bij{x, D) of degree exactly equal to (jj + tj. 

The Lopatinskii conditions are defined as follows. For each x G dX., we denote by 
v{x) the outward unit normal to X at x G dX. We then think of z as the parameteriza- 
tion of the half line x — zv{x) for 2 > 0. Let then C, G S"~^ with C, ■ v{x) = 0. Consider 
then the system of ordinary differential equations 

Ao{x,C + j'^ix)-^)u{z) = in 2; > 
d 

Bq(x,C + ^'^ix)^)u(z) = at 2 = 0. 
* dz 

We assume that the only solution to the above system such that u{z) — t- as 2 — 00 is 
u = 0. This is the Lopatinskii condition for {A, B). We then also say that B covers A. 

The above conditions need to be verified for the specific problems being considered. 
In some situations, the boundary conditions provided by ([1]) (or their linearization) 
generate a cover of A. In other situations, they need to be augmented with additional 
boundary conditions for 5uj as well as for ^7. 

When A is elliptic and B covers A, we say that {A, B) is an elliptic system. 
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2.5 Parametrices and stability estimates. 

Following work in [21 [3l [23] on determined systems, the case of overdetermined elliptic 
systems was treated in |13]. The salient feature of these works is that the operator 
A = {A, B) admits a left-parametrix (a left-regularizer) in the following sense. Let 
(iS, 0) in (E]) be in the space 

7^(p,/) = W^-'^X) X ... X Wl-'^'{X) X wl'^^'^dX) X ... X wl'^^^'^dX), 
for some / > and p > 1 and let us assume that {A, B) is a bounded operator from 

V G v{j), I) = X ... X wj,+''+^'{x) 

to {A,B)v = (iS, 0) G 7l{p,l). Such is the case when the coefficients of A and B are 
sufficiently regular. More precisely, with I sufficiently large so that p{l — Si) > n for 
all 1 < z < 2J (to simplify; see [H] for slight generalizations), we assume that Aij is 
a sum of homogeneous operators of degree Si + tj — n for < k < Si + tj and that 
the coefficients of these operators are of class Wp~''^'{X). Moreover, assuming / large 
enough so that p{l — ag) > n as well for 1 < g < Q, we assume that Bgj is a sum of 
homogeneous operators of degree ag + tj — k ior < k < aq + tj and that the coefficients 

of these operators are of class W''~'^'''p{dX). Here Wp{X) is the standard Sobolev space 
of functionsv with s derivatives that are p-integrable in X with standard generalization 
for s not an integer |T]. 

The main result in [H] is the existence of a bounded operator R from TZ{p, I) to 
©(p, I) such that 

RA = I -T, (8) 

where / is the indentity operator and T is compact in V{p,l). When 1 is not in the 
spectrum of T so that / — T is invertible, then A is invertible with bounded inverse 
{I — T)~^R. However, 1 could very well be in the spectrum of T, in which case dimKerA 
is finite but positive. 

Moreover, we have the following stability estimate 

J+M 2 J Q 

j=l i=l 1=1 p ^ ' tj>0 

for some constants C > and C2 > 0. 

The presence of C2 > indicates the possibility that A may not be invertible. The 
presence of finite dimensional kernels is a serious difficulty in the analysis of the nonlinear 
problem (|3]) because such a dimension is not stable with respect to perturbations. All 
we know is that for Ai sufficiently small, then dimKer(yl + ^i) < dimKerA; see, e.g., 
[28]. 

Whether we can choose C2 = above, i.e., whether A is invertible, depends on 
lower-order terms that are not captured by the principal part {Ao,Bo). Their analysis 
can prove quite complicated in practical settings and we do not follow that route here. 
Instead, our aim is to modify A so that a unique continuation principle may be applied. 
In section [3], we augment the properly modified system {A, B)v = (5, 0) with additional 
boundary conditions, which in some cases allow us to obtain injectivity results. 
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Note that the parametrix R is clearly not unique. It is theoretically constructive, 
as can be seen by following the proof in [3l HI]. However, its practical, for instance 
numerical, implementation is not straightforward. The modified, higher-order, systems 
proposed later in the section offer a more direct numerical inversion procedure. 

2.6 Example of power-density measurements 

To illustrate the theoretical result of this paper, we consider the example of the re- 
construction of a scalar coefficients from the so-called power density measurements. 
Consider the scalar elliptic equation 

£(7, Uj) := V ■ -fVuj = in X, Uj = j) on dX, I < J < J. (10) 

Here, X is an open domain in M" for > 2 with smooth boundary dX. The objective is 
to reconstruct the scalar coefficient 7, uniformly bounded above and below by positive 
constants, from knowledge of the power densities 

Hj{x) = M{n,Uj) ■.= -i{x)\Vuj\^, xeX, 1<J<J, (11) 

where Uj is the solution to (fTOj) . 

This problem and some variations have received significant theoretical and numerical 
analyses in recent years; see, e.g., [5l [TJ HI [131 IIHl IMl EQl [311 [36]; generalizations for 
anisotropic coefficients 7 can be found in [H7]. Explicit reconstruction procedures exist 
when the number of internal functionals J is sufficiently large; see [HI [T^l [32]- The case 
J = 1, which does not correspond to an elliptic system, was analyzed in [7]. The main 
features of this analysis are recalled in section 12.6.11 For intermediate values of J, the 
above hybrid inverse problem may not have an explicit reconstruction but may still be 
modeled by a redundant elliptic system. Such a problem was also analyzed in [301 [H] • 
The conditions of ellipticity of the system [A, B) are described in detail in section 12.6.21 
A modified system is presented in I2.6.3[ whereas optimal stability estimates of the form 
(|9|) are presented in section 12.6.41 for the power density measurement problem. 

2.6.1 The 0— Laplacian when J = 1 

When J = 1, the 2x2 system of nonlinear partial differential equations is formally 
determined with two unknown coefficients (7,^1). The elimination of 7 from such a 
system is in fact straightforward and we obtain the equation for m := Ui (with H := Hi) 
given by 

V ■ T^^Vu = in X, u = f on dX. (12) 

The above equation may be transformed as 

f^ii 

(/-2Vm^Vm) : V\ + Vlni/- Vm = in X, u = f and ^=J on dX. (13) 

Here Vm = and we introduced Cauchy data on dX anticipating the fact that (1131) is 
a quasilinear strictly hyperbolic equation, at least provided that Vti is defined. Indeed, 
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we observe that the operator (J — 2Vm ® Vm) : is hyperbohc with respect to the 
(unknown) direction Vm. 

The above problem is analyzed in [7]. The two salient features of that analysis are: 
(i) unique reconstructions of u, and hence 7, are guaranteed only on part of the domain 
X] see [7]; and (ii) the stability estimates are sub-elliptic: first-order derivatives of u are 
controlled by the gradient of H rather than second-order derivatives as would be the case 
if (/ — 2Vm ® Vm) : was replaced by an elliptic operator. For 7, this translates into 
an inequality of the following form. Let H and H be two measurements corresponding 
to the pairs (^,7) and u,^), respectively. Assume that the Cauchy data of u and it 
agree on dX. Then we find that on an appropriate (see [^) subdomain O (Z X, we have 
the following stability estimate: 

II7 - 7||l^(o) < C\\VH - VHWl^o)- (14) 

As we shall see below, this estimate is sub-optimal (with a loss of one derivative) when 
compared to elliptic estimates of the form Q. It is however optimal for (principally 
normal) operators of principal type |26j . 



2.6.2 Linearization and ellipticity 



The linearization of the above problem with J = 1 is also clearly a hyperbolic equation. 
We now want to show that redundancy in the data (J > 2) allows us to render the 
system elliptic under appropriate conditions. We consider two ways to obtain elliptic 
systems of equations. 

We first linearize the coupled system f|TO|) - f|TT]) about solutions {'y,Uj) and obtain 



V ■ + V ■ -fV6uj 

6uj 





6H, 







in X 
in X 
on dX. 



(15) 



We define 



Vuj and assume that 



Fj I is bounded from below by a positive constant 



uniformly. Such an assumption is valid for and appropriate open set of boundary condi- 
tions fj] see, e.g., [H [H [HI [15] for details of constructions based on complex geometric 
optics solutions or unique continuation principles, which we do not reproduce here. 

Let Aj be the operator applied to v = (7, {uj}) in the above system. Its principal 
part Vj has for (principal) symbol pj a 2 J x ( J + 1) matrix given by 

/ I f?J2 Ors,T?, . n \ 



\Fi\' 27Fi-ze 



\Fj\' 
Fj-ii 











27FJ ■ li 



(16) 



When J = 1, the determinant of pj is given by 7|Fip(2(Fi • E)"^ — |^P), which is 
hyperbolic with respect to Fi as seen in section 12.6.11 above. The above system is in 
Douglis-Nirenberg form for S2k = 1, S2k+i = 0, ti = 0, tj = 1 for j > 2. 
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When J > 1, we observe that sub determinants with the largest number of powers 
of are of the form 

7^|er(-^-^)|F,|\-(x,0, := 2(F,(x) ■ 0^ - 1^1'- 

We thus obtain that pj(x,^) is injective if the quadratic forms qj{x,$,) = for all 
l<j<J imply that ^ = 0. 

Definition 2.1 Define the quadratic forms and operators 

q,{x, = 2{F^ ■ 0' - ler, Pj{x, = A - 2F, ® : V ® V. (17) 

Here Fj{x) are vector fields of unit vectors defined on X. We say that the family {qj} 
or {Pj} is elliptic at x if 

qj{x,C,) = for all 1 < j < J implies ^ = 0. (18) 

We say that such families are elliptic in X if they are elliptic in at all points x G X. 

We can then prove the 

Lemma 2.2 We assume that Fj := Vuj is such that \Fj\ is hounded from below by a 
positive constant on X for all 1 < j < J . 

The operator Aj defined in ( fT5l) with principal symbol given in ( fT6ll is elliptic in X 
if and only if the above family of quadratic forms {qj} is elliptic in X . 

Proof. We have already seen the sufficiency of the condition. Let us prove its 
necessity and assume that pj(x, ^) is maximal rank for x E X and 7^ 0. This means 
that one determinant of ( J + 1) x ( J + 1) sub-matrices of pj is non-vanishing. Each 
column of pj beyond the first one has two non- vanishing entries. For all 1 < j < J except 
for one entry jo? then either — 7(x)|.^p or 2'-)Fj ■ i^ appears as a multiplicative factor in 
the determinant of the sub-matrix. Since 7(x)|.^p never vanishes, we may discard the 
determinant involving '2'^Fj ■ i^. We thus obtain that if one determinant of a sub-matrix 
does not vanish, then that determinant may be chosen as 'y"^\$,\'^^'^~^'^\Fj^^\'^qjg{x, Since 
by assumption 'y^\^\'^^'^~^^\Fjg\'^ is bounded away from 0, we observe that the injectivity 
of pj implies that (at least) one of the quadratic forms qj{x,C,) does not vanish. And 
thus, Aj is elliptic implies that {qj} is elliptic in the sense of definition 12.11 □ 

The ellipticity of Aj is thus a consequence of the fact that the null cones of quadratic 
forms intersect only at 0. We have the following properties: 

Proposition 2.3 (i) Let x G and assume that Fi{x) and F2{x) are neither parallel 
nor orthogonal. Then the corresponding {qi,q2) in- (pTl) form an elliptic family at x. 

(a) In dimension n > 2, let Fi[x) and F2{x) be two different directions and define 
i^3(x) = aFi{x) + PF2{x) with af5 ^ such that \F^{x)\ = 1. Then the corresponding 
{Qi, 92,13) in f|T7|) form an elliptic family at x. 

(Hi) In dimension n > 3, a family (91,52) is never elliptic at a given point x inde- 
pendent of the choice of Fi,F2. 
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Proof, (i) In dimension n = 2, it is clear that the null cones (where qj vanishes, two 
lines of vectors in M^) coincide if and only if Fi and F2 are either parallel or orthogonal. 

(iii) The result is obvious if Fi = ±^2. Assume otherwise and let -F3 be a unit 
vector orthogonal to Fi and F2. Assume Fi ■ F2 > for otherwise change the sign 
of F2 (which does not modify the quadratic forms qj). C, belongs to the intersection 
of the null cones {qj{x,0 = 0} for j = 1,2 if 2(Fi ■ 0^ = 2(^2 ■ 0^ = I^P- Define 
^ = F1 + F2 + AF3. The first constraint is satisfied when 2(1 + A -^2)^ = \F^ + F^\^ + X^, 
i.e., if A = ±(2Fi ■ F2 + 2(A ■ ^2)2)! But then it is clear that this ^ 7^ also belongs 
to the second cone so that {qi,q2) is not elliptic. 

(ii) Let us assume that F3 = aFi + /3F2 with a(3 7^ 0. Then belongs to the three 
null cones if 

Expanding the first constraint, we get 

a\F, ■ 0' + I3\F2 ■ 0' + 2a/3F, ■ ^4 ■ ^ = ^lel'(«' + + 2af3F, ■ F2). 

The last two constraints imply that 2Fi ■ ^^2 ■ ^ = ^ICP with e = ±1. Combined with 
the latter, they also imply that 

a/3e = a/3Fi ■ F2. 

Since a/3 7^ 0, this can only occur if Fi = ±F2, which is a contradiction. □ 

From the practical point of view, this result says that if Fi = Vui and F2 = Vu2 
are not parallel, then the internal functionals Hj for Ui, U2, and Ui + U2 (with F3 = 
V(mi + U2) = Vmi + VM2) generate three quadratic forms qj, I < j < 3 that form 
a elliptic family. Such internal functionals are obtained by choosing three boundary 
conditions of the form /i, /2, and /s = /i + /2- This result holds for all n > 2. 

Sufficient conditions for ellipticity. We have thus obtained the following result. In 
dimension n = 2, J > 2 is necessary for Aj to be elliptic. Moreover, J = 2 is sufficient 
when n = 2 if Vui and Vn2 are nowhere parallel or orthogonal. In dimension n > 3, 
J > 3 is necessary for to be elliptic. Moreover, J = 3 is sufficient for Aj to be 
elliptic in all dimensions n > 2 by choosing as boundary conditions, e.g., (/i, /2, /i + /2) 
provided that Vmi and Vm2 are nowhere parallel. 

Boundary conditions and Lopatinskii condition. In order to obtain an optimal 
theory of stability estimates, the system needs to be augmented with boundary condi- 
tions that satisfy the Lopatinskii condition. Dirichlet conditions on 6uj and no condition 
on ^7 satisfy such conditions. Indeed, we need to show that v{z) = (57(2), . . . , 6uj{z)) = 
is the only solution to 

duj{0) = 0, Fj-Nd,S-f + -fdlSuj = 0, \Fj\^S-f + 2-fFj ■ Ndjuj = 0, z > (19) 

with Suj{z) vanishing as 2; —t- 00 for = ^{x) at x G dX and z coordinate along — A^. 
We observe that this is the case if {Fjl"^ — 2{Fj ■ N)"^ = \Fj\'^qj{x, N) 7^ for some j. 
This is the condition for joint ellipticity described in definition 12.11 
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Let us define the spaces X = W^X) x Wl+\X;W) and y = W^(X;M2J) with / 
large enough so that the latter spaces are all algebras; i.e., pi > n. We also define 

= Wp ''{X■,M^^) for the traces of 6uj on dX, which vanish by construction. Then 
we observe that {A, B) in (IT^ maps V{p, I) = X to n{p, I) = y X y^. Moreover, the 
coefficients 7 and Uj appearing in the definition of A belong to W^iX) and W^^+H^; ^'^)) 
respectively, by assumption for 7 and elliptic regularity for Uj solution of the elliptic 
problem (fTO|) . 

2.6.3 Elimination and ellipticity 

The above system involves J + 1 unknowns. Strategies to lower the dimension of the 
system include: (i) eliminating 7 as we did in obtaining f[T^ : or (ii) eliminating all Uj. 
The second strategy follows from the observation 

6uj = L~\V -S-fVuj), L^ = -V-7V, 

with defined by solving with vanishing Dirichlet conditions. 

The result is a redundant system of the form PjS'j = 5Hj where Pj properly defined 
is a pseudo-differential operator with principal symbol given by \Fj\'^qj{x,^). The re- 
dundant system for ^7 is therefore elliptic under the same conditions as those for (1161) 
above. The main difficulty is that Pj is no longer local (no longer a (system of) partial 
differential equation). It is therefore not clear how one may approach the question of 
uniqueness for such a system. The unique continuation principles presented in section 
[3] do not apply. See [31] for an analysis of such a method. 

The first strategy based on the elimination of ^7 preserves the differential structure 
of the original system since ^7 appears undifferentiated in the second equation in f[T^ . 
We find that 

V ■ 7( - V5u, + 2F,F, ■ V5u,) = V ■ ^F,. (20) 

The symbol for such an equation is then given by 'yqj{x,^) so that the above equation 
is not elliptic, as we already know. However, the elimination of ^7 also provides the 
constraint 

1 6H 1 SH 

27-— F,- ■ V6u, - ^ = 27^i^/c ■ V6u, - I < j < k < J. (21) 

I jl \ j\ \ k\ 

It turns out that the combination of fl20|) with f l2T|) makes the redundant system for the 
{uj} elliptic provided qj{x,^) = for all j implies that = as above. To see this, 
assume that 

qj{x, ^)5uj = 0, \Fk\'^Fj ■ i5uj = \Fj\'^Fk ■ ^5uk. 

Let ^ 7^ 0. Then, not all qj{x,^) vanish. Assume that qi{x,^) 7^ 0. Then 6ui = so 
that Fj ■ ^6uj = for all j > 2. However, Fj ■ ^ and qj{x,C,) cannot vanish at the same 
time so that Suj = 0. This shows the injectivity of the symbol of the redundant system, 
which is easily found to be in Douglis-Nirenberg form. 

Moreover, unlike the redundant system of strategy (ii), the above system fl20|) - (l2Tl) 
is differential. It may therefore be augmented with boundary conditions that satisfy the 
Lopatisnkii conditions. We leave the details to the reader to verify, as we did in ( IT9l) . 
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that such conditions are satisfied for Dirichlet conditions 5uj = on dX under the same 
conditions guaranteeing elhpticity in X. The above system for the {{Suj}) is therefore 
elhptic when the family of quadratic forms {qj} is elhptic. 

2.6.4 Stability estimates 

Both systems of differential equations for (57, {Suj}) and for {{6uj}) after elimination 
of 57 are elliptic with Dirichlet conditions for 6uj under the conditions stated, e.g., 
in Proposition 12.31 We may therefore apply the general theory of elliptic redundant 
systems described above and obtain a parametrix for both systems. Moreover, in both 
cases, we obtain the following stability estimates 

J J J 

<CY,\m,-6H,\\ W{,{X)+C2 \\6Uj-6Uj\\Lp(^x)- 
j=l j=l j=l 

(22) 

In other words, ^7 and VSuj are reconstructed from 6Hj with no loss of derivative, 
unlike the construction in (I14p . However, we are not guaranteed that any of the linear 
systems is indeed invertible, and hence the presence of the last term on the right-hand- 
size in (122|) . Proving the injectivity of the above systems is much more delicate than 
proving their elhpticity. The following section proposes some generic strategies to do 
so. 

2.6.5 Generalization to similar models 

The results presented above generalize to the setting where the internal functionals are 
of the form 

^.(^)=7lVn,f, (23) 

where a > 0. The case a = 1 was treated above. The case a = 2 corresponds to the 
setting of CDII and MREIT SB]. As shown, e.g., in [HB ES], the coupled problem 
is elliptic with J = 1 for a > 2, is degenerate elliptic for a = 2, and is hyperbolic for 
a < 2. The cases a > 2 become elliptic for J chosen sufficiently large. We do not 
consider these extensions further here. 

2.7 Sufficient conditions for ellipticity 

To summarize the above derivation, we have modeled the nonlinear hybrid inverse prob- 
lem as a coupled system of nonlinear partial differential equations (|3]). Its linearization 
is then given by @. These systems have J + M unknowns for 2J > J + M equations. 
When J is large compared to M, we expect the principal part ^0 of -4, in the Douglis- 
Nirenberg sense to be elliptic since a redundant matrix is "more likely" to be full-rank 
than a less elongated matrix. 

The matrix ^0 is full rank if we can prove that, for well chosen boundary conditions 
fj, the internal functionals Hj are sufficiently independent. Such a result is problem- 
dependent. In the setting of power-density measurements, we have obtained that the 
boundary conditions fj were well-chosen if the quadratic form qj{x, ^) were jointly elliptic 
at each point x E X. A sufficient condition to do so is to choose two boundary conditions 



12 



/i and /2 such that Fi = Vui and F2 = Vu2 are nowhere co-hnear (for the quadratic 
forms corresponding to the three boundary conditions /i, /2, and /i + /2 are then jointly 
elhptic). 

Strategies to ensure that, e.g., Vmi and Vu2 are nowhere co-hnear have been pre- 
sented in, e.g., [3 m Uni [E], see also the review [S]. Such strategies are based on the 
use of Complex Geometrical Optics when the latter are available, or on the use on local 
construction and unique continuation properties of the operators in as described in 
[14] . In both settings, it is proved that qualitative properties of solutions to elliptic 
equations, such as for instance the independence of gradients of solutions, hold for an 
open set of boundary conditions fj. 

3 Injectivity results for the linearized problem 

In this section, we consider two methods to obtain injectivity of an elliptic operator 
A augmented with appropriate boundary conditions B. Both are based on replacing 
the redundant system A by its normal, determined, form A* A, and augmenting it 
with appropriate Dirichlet boundary conditions that ensure its injectivity under certain 
assumptions. The first method invokes a Holmgren unique continuation principle while 
the second method is based on unique continuation principles that are consequences of 
Carleman estimates. 

Here, the operator A^ is defined such that {A^)ki = {AikY ='■ A^^., where {AikY is 
the formal adjoint to Aik for the usual inner product (-, ■) on L^(X). 

One difficulty with operators A that are elliptic in the DN sense is that the normal 
operator A'^A need not be elliptic, even in the DN sense. Consider for instance the 2x2 
system in one independent variable defined by = a, A12 = A21 = dx and A22 = d"^, 
which is DN elliptic with si = ti = and S2 = ^2 = 1 when a ^ 1. Defining C = A'^A, 
the principal term in C is Cu = d^, C12 = C21 = d^, C22 = d^, which is independent of a 
and not elliptic. 

One way to ensure that A^A is elliptic when A is elliptic, is to assume that r = 
tj independent of j and Si = independent of i. We then verify that {A'^A)ki is a 
polynomial of degree equal to 2r and that A^A is a strongly elliptic system of size 
{J+M)x{J+M); see, e.g., [IZl Proposition 4.1.16]. The principal part of Det(>l*^)(a;, 
is equal to Det(^Q^o)(^! and is a polynomial of degree 2( J + M)r that is uniformly 
bounded from below for all x G X and ^ G S"~^ by assumption of ellipticity. 

The results in [20] show that an elliptic system in DN form can always be transformed 
into an elliptic overdetermined system of first-order equations. The procedure increases 
the order of the system and differentiates any row involving a term with Si + tj = 0. 
This forces us to impose boundary conditions on the parameters ^7 that may not be 
necessary in the definition of {A, B) above. 

For the rest of the paper, we assume that A has been recast into a form where r = tj 
is independent of j and Sj = is independent of i; for instance with r = 1 or r = 2 as 
described in [20]. We still keep the notation 2 J and J + M for the size of the system A 
so that with this notation, A'^A is a system of size (J + M) x ( J + M). 

Let us augment A'^A with Dirichlet boundary conditions 

d 

{jT-Yvj = (j)gj ondX, < g < r - 1, I < 3 < J + M. (24) 
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We recast the above constraints as Vv = (j) on dX. It is proved in [21 p. 43-44] that such 
boundary conditions cover A'^A, i.e., that the Lopatinskii conditions are satisfied. We 
thus consider the problem 

A*Av = A^S inX, Vv = (p ondX. (25) 

Since := {A^A,V) is eUiptic, the above system admits a left parametrix G such that 
GN = I — T with T compact in V^p, I). Moreover, we have the stability estimate 

J+M J+M J+M 

(26) 

We verify that 

2J 2J 

i=l i=l 

so that ( |26ll may also be seen as an analog of ([9]). 

Our objective is to find sufficient conditions under which the above system is injec- 
tive, and hence invertible, so that C2 = in the above estimates. We start with the 
following simple lemma: 

Lemma 3.1 Let us assume that v is a solution of 

A^Av = mX, Vv = on dX. (27) 
Then v is a solution of 

Av = m X, Vv = on dX. (28) 

Note that the Dirichlet conditions are now redundant for ( l28l) . For example, consider 
^ = A + c(x) a scalar operator. For some choices of c(x) (for instance c{x) = A an 
eigenvalue of —A on X with Dirichlet conditions), A is not invertible. However, fl28|) 
corresponds to (A + c{x))v = with both v = and d^v = on dX. It is then 
known that the unique solution to the above constraints is f = and hence {A^A, T)) is 
injective. 

Proof. We observe that the differential operator {A^)ki = {AikY ='■ Ajj. is of the 
same order r as Aik- Since d^Vk = for < g < r — 1, we obtain by integrations by 
parts that 

= Y(AlA,v,,Vk) = = Y 11(^^)^11'- 

ijk i i k i 

This implies the result. □ 

3.1 Holmgren unique continuation and generic results 

The above result shows that injectivity of ( l25l) is a consequence of the injectivity of 
(128!) with redundant boundary conditions. The operator A is therefore augmented with 



14 



more boundary conditions than is necessary to render it elliptic. This is a possible price 
to pay to guarantee that the solution to the normal system (!25|) is uniquely defined. 

In general, however, even with redundant boundary conditions, it is not entirely 
straightforward to ensure that w = is the only solution to ( 128|1 . In this section, we 
consider the case where A is well-approximated by an operator with analytic coefficients. 
We have the following result adapted from |27j . 

Proposition 3.2 LetAAix,D) be a2Jx[J+M) elliptic system of differential equations 
with analytic coefficients in X. We assume that J + M < 2J. Moreover, let us assume 
that Aa{x,C,) is of full-rank J + M for every x & X and ^ G Assume that Aav = 

in X . Then v is analytic in X. 

Proof. Since Aav = 0, we also have that 

^\A\Aa){A\Aa)v = Diag(det(^*^^A), • • • , det(^^^A))f = 0, 

where '^"{A^a-^a) is the matrix of co-factors of A\Aa- Since Aa is of full rank for every 
^ e then P = det(^^^yi) is an elliptic operator such that CharP = 0, where the 

characteristic set of P is {(x, ^) G X x Pm{x,^) = 0} and Pm is the principal part 

of P. We deduce from [271 Theorem 6.1] that WFA{vj) C CharP U WFA{Pvj) = so 
that each vj is analytic in X. □ 

From this, we deduce the following unique continuation result: 

Theorem 3.3 (Holmgren) Let Aa be as in Proposition \3.B and let us assume that 
Aav = in X. Then we have: 

(i) Assume that v = in an open set Q CC X . Then v = in X . 

(a) Assume that Vv = Q on an open set S of dX. Then f = in X . 

Proof. Let us first prove (i). We know from Proposition 13.21 that the functions Vj are 
analytic. Since they vanish on an open set, they have to vanish everywhere. 

Now (ii) is a standard consequence of (i), which we write in detail for systems. Let 
Xq G S and V a sufficiently small open ball around Xq where the coefficients of Aa 
are analytic and where Vv = on E fl C S. Since Aa is injective, we deduce that 
(dnY^Vj = on S n y as well, as any higher-order derivative in fact. Let us extend v by 
on V. Then we verify that Aav = in ^ (all that needed verification was that the 
equation was satisfied point-wise on T, (1 V) . But now (i) implies that = on ^ U X. 
□ 

The above result shows that Aa = AaA is injective as well. Since Ra^a = I — Ta 
for some left-parametrix Ra and compact operator T^, we deduce that 1 is not in 
the spectrum of Ta so that A']^ = (J — Ta)~^Ra- As a consequence, any operator A 
sufficiently close to Aa is also invertible with a bounded inverse. We summarize this 
result as: 

Corollary 3.4 Let A = Aa + Ai with Ai sufficiently small in operator norm from 
Yli^pi^) 'to Y\jWp^^{^)- Then A augmented with Vv = is injective and A* A 
augmented with the same boundary conditions is invertible with bounded inverse. In 
other words, ( j26l) holds with C2 = 0. 
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The invertibility of A^A is therefore generic, i.e., holds for an open and dense set of 
coefficients in the definition of the elhptic operator A^A; see [15]. Note, however, that 
the size of Ai for which A is invertible depends on Aa- 

Another corollary states that the invertibility of A'^A is guaranteed when X is a 
sufficiently small domain. Indeed, let us assume that all coefficients of A are sufficiently 
smooth and that G X. Let Aoo{D) = A{0, D) be the operator with coefficients frozen 
at a; = 0. Then ^qq^qo is invertible as an application of Theorem 13.31 The operator 
A{x, D) — A{0, D) is then small in operator norm provided that X is sufficiently small. 
As in the derivation of the preceding corollary, we thus obtain that A is injective and 
that A^A is invertible when augmented with Dirichlet conditions Vv = 0. We summarize 
this result as: 

Corollary 3.5 Let us assume that A{0,D) is an elliptic operator. Then for X suffi- 
ciently small, we have that A is injective and that A^A is invertible when augmented 
with Dirichlet conditions Vv = 0. As a consequence, fl2B]l holds with C2 = 0. 

As in the previous corollary, the size of the domain X depends on the bound for 
{A^^qqAqq)^^ ■ The Holmgren result is not quantitative and does not provide any in- 
formation regarding that bound. 

3.2 Unique continuation principle 

When the domain X is large or when the operator A is not sufficiently close to an 
operator Aa with analytic coefficients, then unique continuation results must rely on an 
other principle than analyticity. 

Unique continuation in the absence of analyticity in the coefficients is not always 
guaranteed, even for scalar elliptic operators, although the construction of counter- 
examples is not straightforward. For general references on unique continuation principles 
(UCP) for mostly scalar, not necessarily elliptic, operators, as well as counter-examples, 
see [281 Ell HOl [50] and their references. 

In this section, we revisit Calderon's [TTJUH] result of uniqueness from Cauchy bound- 
ary data closely following the presentation in [IQ] and extend it to redundant systems. 
Our objective is to adapt [1U| Theorem 5] to specific settings of redundant systems. We 
first present local uniqueness results, whose proofs are postponed to the appendix, and 
then use classical arguments to extend them to global uniqueness results. 

3.2.1 Local uniqueness result 

Let a; be a point in W^^^ and iV be a unit vector equal to (0, ... , 0, 1) in an appropriate 
system of coordinates. We want to address the local uniqueness of the Cauchy problem. 
Assume that Lq ioi 1 < q < Q are differential operators of order m and that 

LgU = in y n {xn+i > 0}, dl^^u = on n = 0}, (29) 

for 1 < q < Q and 1 < j < m — 1, where ^ is a neighborhood of 0. We assume that 
is non characteristic (at x) for all Lg, i.e., Pq{x, N) ^ with pq the principal symbol of 
Lq. This and (129|) imply that all derivatives of u vanish on {xn+i = 0} and that u can 
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be extended by on n {xn+i < 0}. The uniqueness problem may therefore be recast 
as: if u satisfies 

LqU = in V ioT 1 < q < Q and m = in {xn+i < 0}, (30) 

then M = in a full neighborhood of 0. 

When Q = 1, it is known that L = Li needs to satisfy several restrictive assumptions 
in order for the result to hold; see [ID]. The main advantage of the redundancy in the 
above system is that such assumptions need to be valid only locally (in the Fourier 
variable C,) for each operator, and globally collectively. 

Changing notation, we define t = Xn+i and still call x = [xi, . . . ,Xn)- Then 
Pq = Pq{x,t,^,T) is the principal symbol of Lg for this choice of coordinates. Suffi- 
cient conditions for the uniqueness to the above Cauchy problem involve the properties 
of the roots r of the above polynomials r (-> Pq{-, t) as a function of ^. In the setting of 
redundant measurements, these conditions may hold for different values of q for different 
values of ^. This justifies the following definition of our assumptions. 

We assume the existence of a finite covering {Qu} of the unit sphere (corre- 
sponding to 1^1 = 1) such that the following holds. For each z/, there exists q = q{u) 
and e > such that for each (x, t) close to and each ^ G fi^, we have: 

(i) Pq{x,t,^,T) has at most simple real roots r and at most double complex roots 
(a) distinct roots Ti and T2 satisfy |ti — > s > 
{Hi) non-real roots r satisfy \'^t\ > e. 

Then we have the following result: 

Theorem 3.6 [Calderdn's result for redundant systems.] Assume that N is non char- 
acteristic for the operators Lg at the origin and that for a finite covering {fi^} of the 
unit sphere, above are satisfied. Then (l30l) implies that u = in a full 

neighborhood ofO. 

The proof, which closely follows that of Theorem 5 in [ID], is presented in the appendix. 

In the above theorem, u is scalar. The above proof extends to vector- valued functions 
u when the system for u is diagonally dominant. For a determined system given by a 
matrix Lij of operators of order m > 1, this means that the operators are of order 
m and satisfy the hypotheses of the above theorem and the operators for i ^ j are 
of order at most m — 1. This generalizes to the setting of redundant systems LijUj = Si 
where for each i, only one operator Lij is of order m and for each j, the operators L^j 
of order m collectively satisfy the hypotheses of the above theorem; see Theorem IA.2l in 
the appendix, which we do not reproduce here. 

The above theorem may not apply directly in applications. However, it serves as 
a component to obtain more general results. We consider one such result that finds 
applications in the framework of power density measurements. 

We consider a setting where the leading term in the system may not be diagonal 
but rather upper-triangular. Unique continuation properties may still be valid provided 
that (a sufficiently large number of) the diagonal operators are elliptic. However, the 
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corresponding complex roots may no longer be double. Instead, we need the stronger 
condition for some 5 > 0: 



(iv) Pq{x,t,^,T) has at most simple roots r that satisfy |$5r| > e. 

Then we have the following result. 
Theorem 3.7 Consider the redundant system of equations 




= zn y n {xn+i > 0}, di^^Uk = onVn {xn+i = 0}, (31) 



for 1 < j < m — 1 and 1 < A; < 2, with the following assumptions. The operators Li and 
Lq are (vector-valued) Qixl operators of order m, where Li satisfies the hypotheses of 
Theorem \3. 6\ with Q = Qi. The operators L2 and L3 are Q2 x 1 operators of order m and 
at most m — 1, respectively. Moreover, L2 satisfies the ellipticity hypothesis (i)-(ii)-(iv) 
with Q = Q2. Then {ui,U2) = in a full neighborhood ofO. 
The same result extends to systems of the form 

(Lij) j < = in Vn{xn+i>0} , di_^_^Uk = on Vn{xn+i = 0} , (32) 

for 1 < k < R, where Lij is (a vector-valued operator) of order m — 1 for i > j , La is 
(a vector-valued operator) of order m that satisfies the hypotheses of Theorem \3.b\ with 
an appropriate value of Q when all Li^ are of order m — 1 for k ^ i and the ellipticity 
hypothesis (i)-(ii)-(iv) with an appropriate value of Q when at least one operator Ln. for 
k < i is of order m. 

The proof of this theorem is also given in the appendix. Note that (i)-(ii)-(iv), as 
opposed to (i)-(ii)-(iii), may be seen as an ellipticity condition since the symbol Pq(^^) 
does not vanish on Q^,. It is known that Carleman estimates are sharper in such a 
setting; see, e.g., [321 US]- 

Remark 3.8 Let Q be an invertible Qi x Qi matrix and define the linear operators 
Li = QLi and Lq = QLq. Then the conclusions of Theorem 13.71 are equivalent to 
the those when Li and Lq are replaced by Li and Lq. In other words, the results of 
Theorem 13.61 hold if Lq is replaced by an equivalent linear combination. The condition 
that be not characteristic for Lq may then be replaced by the condition that it be 
not characteristic for the linear combinations Lq = Qq^qiLqi. 

Remark 3.9 If Lq is an operator of order m — 1 in the latter theorem so that the 
above system is diagonally dominant, then both Li and L2 are required to satisfy the 
less constraining hypotheses of Theorem 13. 6j see Theorem IA.2I and the discussion below 
Theorem 13.61 

3.2.2 Global uniqueness result 

The above local results extend to global unique continuation results such as [271 Sections 
8.5-8.6]: 
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Theorem 3.10 Let A be a (redundant) system of operators of order m, let i?^ C 
{{x,N) G X X be the set of points where the hypotheses of either Theorem \3.6[ 

or Theorem 3.1, (and possibly their extensions in Remark \3. 9\) are satisfied and let us 



define = X x Let 

Au = in X, dlu = on dX, < j < m. — 1. 

Let iV(supp(M)) be the subset of (x,^) G X x Ei""^^ composed of the closure of the normal 
set of the support of {uj}; see f2^ . Heuristically, that set may be seen as the points 
{x,^) where x belongs to the boundary o/supp(u) and ^ is a normal (inward or outward) 
unit vector to that boundary. 

Then 7V(supp(n)) C S_4. When the latter set is empty, then u = in X. 

The proof follows from the local results noting that the boundary to the support of u 
cannot occur at a point and for a direction where a local UCP principle holds; see 
Sections 8.5-8.6] for the definitions and derivations. 

Remark 3.11 Let ^ be a p x g system with p > q and let us assume that Au = in 
the above theorem is replaced by 

<? 

\{Au)k\ix)<C El^^^'K^)' (33) 

|a|<m-l 1=1 

for a constant C independent of x G X and 1 < k < p. Let, as above, E_4 be the set of 
points where the hypotheses of Theorem 13.61 do not hold (with obvious modifications for 
Theorems 13. 7p . Then the conclusions of the above theorem, namely that iV(supp(n)) C 
E_4, still hold; see, e.g., [2H] and the proofs in the appendix. 



3.3 Application to power-density measurements 

We now apply the unique continuation results of sections 13.11 and 13.21 to the setting of 
power density measurements described in section 12.61 

3.3.1 Second-order linear systems of equations 

Both unique continuation results require that the system A be elliptic in the regular 
sense, i.e., with r = tj independent of j and Sj = for all i. The systems f|T5|) and 
(l20l) - fl?Tl) are not in this form. We propose two modifications of the latter systems for 
which uniqueness results can be proven. 



System after elimination. Let us first consider the system fl2U]) - fl2T]) . The latter 
constraints are first-order. However, they easily become second-order by differentiation. 
We thus consider the system 

V.7(-V<5n,+2F,F,-V5n,) = V ■ ^F,. 



V(27-^F,.V5«,-27-^F,.V5«.) = V(^ - ^) 1 < j < k < J, 



(34) 
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in X with boundary conditions 5uj = d^Suj = on dX. The above system is clearly el- 
liptic under the same hypotheses as (!20l) - (l2Tl) . Moreover, it satisfies the above ellipticity 
condition with r = 2. 



System in triangular form. Let us come back to (ITS!) . We apply the operator 
Kj = 2'y\Fj\~^Fj ■ V to the first line and = V ■ 7V|-Fj|^^ to the second line to obtain 
after subtraction 



V ■ + V • -fV6uj = 
L^6'j — KjV ■ 6'yFj + [L^Kj, KjL^]6uj = 



SH, (35) 



1^. 



|2 



Since is second order and Kj is first-order, the commutator [L^Kj, KjL^] is a second- 
order differential operator. We thus again obtain a system that is elliptic when (fT5|) is 
(we leave this proof to the reader; see also (l38l) below) and that is in the appropriate 
form with r = 2. 



3.3.2 Holmgren unique continuation result 

The results of section 13.11 apply to the power density measurement systems fl34p and 
(l35l) . Let A denote the linear operator in one of the aforementioned systems. Then A^A 
augmented with Dirichlet conditions is injective provided that 7 and uj are sufficiently 
close to analytic coefficients or provided that the problem is posed on a sufficiently small 
domain, and provided of course that the quadratic forms {qj} are collectively elhptic on 
X in the sense of Definition 12.11 

The main difference between (l34l) and (!35|) pertains to the boundary conditions. In 
the case of (l35l) . vanishing Dirichlet and Neumann boundary conditions are imposed 
on S'y and 6uj. This means that the value of (7, {uj}) and of its (normal) derivatives 
is known on dX. In the case of only the values of {{uj}) and of its (normal) 

derivatives needs to be known. 

In fact, the boundary conditions for ^7 at x G dX may be deduced from those for 
6ui, say, provided that z/(a;) is not characteristic for Pi. Indeed, we observe that 

This provides boundary conditions for ^7 when 6ui and dy5ui are known on dX. The 
derivation of 0^5'') is as follows. Using the equation 

V ■ 57VMJ + V ■ 7V5mj = 0, 

and replacing ^7 above using the expression of we obtain after straightforward 
eliminations that 

where /i is a known expression involving d^Ui and tangential derivatives of Ui of degree 
up to 2 at X, which are known since Ui is known on dX. Since v is not characteristic for 
Pi, this provides an expression for c^^Mi, and hence for (9i/(57 using the above expression 
for 5 Hi. 

So knowledge of Cauchy data for 57 is in fact a consequence of knowledge of Cauchy 
data for the solutions 5uj. Thus, the major requirement in solving A^Av = A'^S is that 
we impose that dy5uj = 0, which implies that d^Uj is known on dX. 
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3.3.3 Calderon unique continuation result 

The symbol of the operator (ITS!) is elliptic in the Douglis-Nirenberg sense. Recall that 
the results of unique continuation principle (UCP) of the preceding sections apply to 
operators that are elliptic in a classical sense. We thus consider the modification of the 
system We were not able to obtain a satisfactory global UCP using Carleman-type 
estimates for the system flMl) . 

UCP for such the modified system fl35l) is based on proving a collective UCP for 
quadratic forms, which we now address in detail. 

Lemma 3.12 Let Pj{x,D) be second-order operators with principal symbols given by 
qj{x,^) = 2{Fj ■ 0^ - |<eP, where Fj{x) E S""^ for 1<J<J. Let {x,N) E X x 
For C E M" such that ^' ■ N = 0, we define 

O = 2(F, ■ - (1 - 2(F, ■ iVf ) lef. 

Let us assume that the quadratic forms are independent in the following UCP sense: 

gj;Ar(x,^') = for all 1 < j < J s.t. qj{x,N) ^ implies that ^ = 0. (36) 

Then {Pj} collectively satisfies a UCP for {x,N) in the sense that we can find a 
partition {Q^} such that (i)-(ii)-(iii) holds. 

Proof. The proof in dimension n = 2 is trivial: Pj satisfies UCP except for four 
directions on the null cone, i.e., for ^ such that qj{x,^) = 0. Since ( l36l) implies that 
is not in the intersection of the null cones, then UCP holds for one of the operators Pj 
and we can choose Qi = §^ in the partition. 

Let (x, N) EX X and n > 3. Let ^' E §"~^ such that ^' ■ = 0. Then 

qj{x,^'+rN) = 2{Fji^'+TN)f-\^'+rN\^ = {2{Fj-Nf-iy+AFyNFj-^'T^-2{Fj-Cf-l 

This quadratic polynomial (because qj{x,N) = 2{Fj ■ N)"^ — 1 7^ 0) with real- valued 
coefficients has a double root when the determinant 

a;. = -1 + 2((F, ■ NY + (F, ■ O') = g,;7v(x,0 = 0. (37) 

For C,' away from such points, then the UCP condition holds for Pj since the roots 
cannot be simple and complex roots (which come in conjugate pairs) cannot be close to 
each-other. By assumption, A^- = for all j is not possible and by continuity, maxj \A'j\ 
is bounded away from so that different roots r are separated by some e > (as well 
as the imaginary part of such roots since they come in conjugate pairs). This generates 
a finite partition {^2,^}. On each patch, one Pj satisfies the UCP properties (i)-(ii)-(iii). 
□ 

In dimension n = 2, ellipticity of (Pi,P2) is equivalent to UCP of (Pi,P2) at each 
(x, A^). In dimension n > 3, collective ellipticity and collective UCP are different notions. 
As soon as there is one j such that 2(P,- ■ N)"^ > 1, the UCP is satisfied at (x, ^) while 
ellipticity does not necessarily holds. However, we have the result: 

Lemma 3.13 In dimension n = 2, collective UCP and collective ellipticity of operators 
of the form Pj above is equivalent. In dimension n = 3, collective ellipticity of {Pj} 
such that qj{N) 7^ implies collective UCP at {x,N). 
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Proof. The case = 2 is handled as in the proof of Lemma I3.12[ Consider now 
the case n = 3. Non-UCP at (x, A^) imphes from (137|) and the fact that is not 
characteristic for Pj that 

(F, ■ N)' + (F, ■ a' = I- 

If we decompose Fj = Fj ■ NN + Fj ■ + F/' with F/' orthogonal to span(iV, ^'), then 
we find that 

In dimension = 3, we find that F," = ±|Fj|A^ x This implies that N x ^' belongs 
to the null cone of qj, i.e., qj{x, N x C,') = 0. This holds for every j, which implies that 
X ^' = 0, a contradiction. □ 

Definition 3.14 We say that the family of operators Pj{x,D) or of quadratic forms 
qj{x,C,) collectively satisfy a global UCP on X if for every {x,N) G X x we can 

find a family qj{x,^) = Qjkqk{x^C) with Q an invertible matrix such that {qj{x,^)} 
collectively satisfies a UCP at {x,N) in the sense of Lemma \3.1B, 

The above definition involves a linear change of quadratic forms for each (x, A^) that 
allows us to obtain a family of modified forms qj with similar properties to those of 
qj and such that qj{N) ^ 0; see Remark 13.81 This will prove useful in the analysis of 
reconstruction from power density functionals. We will also need the following Lemma 
in dimension n = 3. 

Lemma 3.15 Let Fi and F2 be three-dimensional vector fields on X such that for each 
X e X, rank(Fi,F2) = 2. Let F3 = Fi + F2. We denote by {Pj} and {q^} the cor- 
responding operators and quadratic forms defined in f|T71) for 1 < j < 3. Then {qj} 
satisfies a global UCP on X . 

Proof. Let us fix [x, N) e X x and choose a basis of M.^ such that Fi (x) = Ci 
and span(Fi,F2) = (61,62). The three quadratic forms qj allow us to obtain by linear 
combination any qc corresponding to Fc = C6i + 562 with s = \/\ — It is then not 
difficult to find three values of c such that the corresponding qj for 1 < j < 3 form an 
elliptic family and such that qj{N) 7^ 0. From Lemma r3.13[ we obtain that {qj} satisfies 
a UCP at {x,N). The hypotheses of global UCP in Definition 13. 141 are met. □ 



Elliptic system in triangular form. We now present a modification of the linear 
system for (^7, {Suj}) for which a global UCP result as described in the above definition 
can be obtained. 



We recast fl35l) as 





let. + Sj, (38) 



where l.o.t. means a system of differential operators of order at most 1. Let P and 
P denote the J x 1 columns of the second-order, homogeneous, operators Pj and Pj, 
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respectively. The symbol of Pj is proportional to qj{x,^). The operator A satisfies 
hypothesis (iv) of Theorem 13.71 while the operator P is of order 2. If the operators P 
collectively satisfy the hypotheses of Theorem 13. 6[ then the above system (135|) satisfies 
a UCP, as does the fourth-order operator A^A as constructed in the preceding section: 



Theorem 3.16 Let Av = S in X he the system for v = {S'y, {Suj}) described in 
and augmented with boundary conditions v = dyV = on dX . Assume that the co- 
efficients in A are sufficiently smooth (see proof of Theorem \A.1\) . Let the operators 



Pj above satisfy a UCP on X as described in Definition \3A^ and Lemma \3.12 . Then 



any solution to that system, or to the system A^Av = A^S with the same boundary 
conditions, is unique. 

Proof. The proof is a direct application of Theorem 13.71 since A satisfies hypothesis 
(iv), P is of order 2, and P collectively satisfies a UCP at each (x, A^) in X x □ 

In dimensions n = 2 and n = 3, we have the following sufficient conditions on the 
internal functionals to guaranty injectivity of A^A: 

Corollary 3.17 Assume that (mi,M2) are solutions such that Fi = Vmi and F2 = Vm2 
are such that rank(Fi, F2) = 2 for each x & X . 

In dimension n = 2, assume moreover that Fi ■ F2 ^ for each x & X . Then 
{Pi,P2} satisfies a global UCP property on X and the results of Theorem \3.16[ hold. 

In dimension n = 3, define F3 = Fi + F2. Then {Pi, P2, P3} satisfies a global UCP 
property on X and the results of Theorem \3.16[ hold. 

The proof is immediate using Lemma 13.151 remark 13.81 and the preceding Theorem. 

We have thus exhibited a system of equations A'^Av = A'^S with boundary conditions 
V = dyV = on dX , which admits a unique solution that verifies the stability estimate 
(1221) with C2 = 0. 

As in the preceding section devoted to the Holmgren uniqueness result, this result 
comes at the cost of having to impose boundary conditions of the form v and d^v to 
both V = Suj and v = S'j. As we saw in the preceding section, the boundary conditions 
for ^7 at X G dX may be deduced from those for 6uj. 

Comparing ellipticity and uniqueness criteria for A in (1381) . we observe that ellipticity 
is obtained by imposing conditions 6uj = on dX and ensuring that {qj} is an elliptic 
family of quadratic forms. UCP is obtained by imposing the additional boundary con- 
ditions dySuj = and by ensuring that {qj} satisfy a global UCP on X. In dimension 
n = 2, global ellipticity and global UCP for {qj} are equivalent. In dimension n = 3, 
both are satisfied for the family of three internal functionals described in Corollary 13. 171 



4 Local uniqueness for nonlinear inverse problem 

4.1 Iterative fixed point and reconstruction procedure 

Let us consider a setting where {A, B) is injective or {A^A, T>) is invertible so that the lin- 
earized problem ([5]) is solved by either v = {I —T)~^ RA{S , 0) or by v = {A^A)^^{A^S, 0). 
We now present standard results concerning the nonlinear hybrid inverse problem ([3]). 
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Let vq = (7, {uj}) G A* = V{p, I) the point about which was hnearized to give 
dl]). We recast both equations in ([3]) as 

T{vo + v) = n, (39) 

where "H G 3^ collects the available data {Hj} and 3^ is defined as the appropriate subset 
of 7l{p, I) defined in section [2751 We also define "Ho := -^(^^o)- 

In the procedure of construction of a left-parametrix to A, we impose in addition 
to (15^ the boundary condition B{vq + v) = B{vq) so that B{v) = 0. Note that in the 
setting of power-density measurements, we do not have injectivity of such an operator 
A. Let us assume nonetheless that 1 is not an eigenvalue of the operator T so that 
(/ — T)~^i? is a bounded left-inverse for A. 

In the inversion of {A^A,V), we impose in addition to fl39|) the boundary condition 
^(fo + v)= T>{vo) so that T>{v) = 0. We have obtained in the preceding section several 
unique continuation results that apply to the differential operator A given in flMl) or 

We recast (1391) as 

r{vo)v = H- J^ivo) - {J'ivo + v)- J^ivo) - T'{vo)v). 

We are in a setting where J^'(fo) augmented with appropriate vanishing boundary con- 
ditions admits a left-parametrix, which we denote by (J-'')~^(fo), and which takes the 
expression (J — T)~^R or {A^A)~^A^ with boundary conditions Bv = Q and Vv = 0, 
respectively. 

We then have that 

V = g{v) := {rr\vo)i'H - no) - {T'r'MiHvo + v)- J^ivo) - T'{vo)v). (40) 

Since J^iw) is polynomial in w with smooth coefficients and we have the stability 
estimate ([9]) with C2 = 0, we obtain that 

\\giv)-giw)\\y < C,i\\v\\y + \\w\\y)\\v - W\\y 
Ur)-\vo){Hvo + v)-J'{vo)-r{vo)v)\\y < Cs\\v\\l, 

for \\v\\y sufficiently small since y forms an algebra. 

For \\v\\y < R sufficiently small and \\{J-'')~^{vo){'H — 'Ho)\\y < r] sufficiently small so 
that rj + C^R^ < R and 2CiR < 1, we find that ^ is a contraction from B = {v, \\v\\y < 
R} to itself. This proves that for "H sufficiently close to Hq, then the solution to (139!) is 
unique and may be obtained by the converging algorithm v'''^^ = Q{v''). 

Provided that {J^'{v))^^ is shown to exist and to be uniformly bounded as an operator 
from y to X for v in the vicinity of Vq, then a Newton scheme with faster convergence 
properties may be introduced: 



{F\v, + v^))-^ {n - T{vo + t;')) . (41) 



The analysis of the scheme and its very fast (quadratic) convergence rate are then 
classical. 
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Application to power density measurements. Consider the scheme introduced 
in (135|) or (134|) . which we recast as Av = S. Then J^'{vq) is an operator from X = 
W^{X) X VT^+H^; K"^) to 3^ = Wj^iX; M^^), which to S associates v solution of 

A'Av = A'S in X, v = d^v = on dX. (42) 

We have presented sufficient conditions in the preceding section for the above problem 
to be uniquely solvable and hence for {J^'{vo))~^ to be bounded from y to X . 

We then find that J^'{v) is invertible from y to X for \\v — vq\\x sufficiently small 
so that the nonlinear iteration schemes in fHOj) or f HT]) converge. This is because in the 
setting of power density measurements, elliptic regularity for vq G X shows that 6v E X 
as well provided that 6H e y. More precisely, for 7 e Wp{X), then Uj e Wp'^^{X] R-^) 
by elliptic regularity of ffTU]) |H] and then (^7, {Suj}) E X by elliptic regularity of either 
( 135|) or (!34|) since power density measurements for 7 G lyp(X) and u G iyp"'"^(X) of the 
form 7|VmP are indeed in Wp{X). 

When one of the unique continuation results applies (which may require that / above 
be quite large so that coefficients are sufficiently smooth; see the proof of Theorem lA.ll 
in the appendix), then the nonlinear hybrid inverse problem may be solved locally using 
(HOl) or fHB . Indeed, for (jH]), v'''^^ — belongs to X as soon as {J^'{vq))~^ is bounded 
from y to X and T-L—l-io is sufficiently small in 3^. This provides a local uniqueness result 
and an explicit reconstruction procedure for the power density measurement hybrid 
inverse problem. 

We refer the reader to (TU] for the implementation of a similar scheme that displays 
very fast convergence rates but for which we are unable to prove a unique continuation 
principle. 

4.2 Injectivity results 

The fixed point algorithm of the preceding section provides us with an explicit local re- 
construction procedure for the nonlinear hybrid inverse problem. A similar methodology 
allows us to obtain local uniqueness results that are more general but not constructive. 
Using the same notation as in the preceding section, let us assume that 

F{vo + v) = :F{vq) = H, (43) 

in other words the measurements associated with vq and Vq + v are identical. Injectivity 
of the nonlinear problem J-" locally in the vicinity of Vq then means proving that = 0. 

Since J-" is a polynomial in v (or more generally a differentiable function in v), then 
there is a polynomial function in v and v such that 

J^{vo + v)- T{vo + v) = g{v, v)-{v-v). (44) 

When V — V is sufficiently small, then Q{v,v) satisfies the same properties as J^'{v), 
the Frechet derivative of J-". For instance, the ellipticity and unique continuation prop- 
erties of J^'{v) still hold for Q{v,v) for v — v small in the appropriate topology. As a 
consequence, we obtain that P3|l implies that f = and hence that J-" is injective. How- 
ever, Q{v, v) may still admit a left-inverse for v and v not necessarily close to each-other, 
in which case we also deduce that v = v for v — v not necessarily small. 
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Let us assume more generally that we have two measurements 

J^{vq + v) = n, J^{vo + v)=H. 

Assume that we have been able to prove that -F'(fo) is elliptic and injective. Then, 
Q{v^ v) is elliptic and hence admits a left-inverse (since it is injective), at least for v and 
V sufficiently close to vq. As a consequence, if we assume that the boundary conditions 
that need to be imposed for v and Vq agree on dX, then we obtain the stability estimate 

\\v-v\\x<C\\H-H\\y. (45) 

Errors in the boundary conditions B{v) or 'Div) may be accounted for in a similar 
manner. Using the notation in section 12. 5[ this is 

J+M 2J 

j=i i=i 

Note that the nonlinear problem (1431) is invertible generically. Indeed, for vq analytic, 
the inversion of A!^A augmented with Dirichlet boundary conditions V is guaranteed. 
As a consequence, the above result shows that P3|) may be inverted for {y + vq) in an 
open set including vq. Since analytic coefficients fo in restricted to X are dense in 
the set of sufficiently smooth coefficients on X, we obtain that the inverse problem may 
be inverted generically^ see |45j . 

When vq is not analytic, then we need to find a unique continuation principle based on 
Carleman estimates; see the paragraph on applications to power density measurements 
in section HTTI which extends to the setting of injectivity results and not-necessarily local 
stability estimates of the form (l45l) . 
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A Unique continuation for redundant systems 
A.l Redundant system for a scalar function 

Let x be a point in ]R"+^ and be a unit vector equal to (0, ... , 0, 1) in an appropriate 
system of coordinates. We want to address the uniqueness of the Cauchy problem. 
Assume that Lq for 1 < q < Q are operators of order m and that 

LqU = QmV n {xn+i > 0}, dl^^u = on V n {x„+i = 0}, (47) 

for 1 < q < Q and 1 < j < m — 1, where \^ is a neighborhood of 0. Here, m is a 
scalar function solution of a redundant system of equations. We assume that N is non 
characteristic for all Lq. This implies, as in section [3. 2. that all derivatives of u vanish 



26 



on {xn+i = 0} and that u can be extended by on \^ fl {xn+i < 0}. The uniqueness 
problem may therefore be recast as: if u satisfies 

LgU = in for 1 < g < Q and m = in < 0}, (48) 

then M = in a full neighborhood of 0. 

When Q = 1, it is known that L = Li needs to satisfy some restrictive assumptions in 
order for the result to hold; see [261 ESI HOI ED] . The main advantage of the redundancy 
in the above system is that such assumptions need to be valid only locally (in the 
Fourier domain) for each operator, and globally collectively in a sense defined below. 
Our derivations closely follow the presentation in [30] , to which we refer to additional 
details. 

Changing notation, we define t = Xn+i and still call x = (xi, . . . , x„). We then define 
Pq = Pqi^y ^, t) as the principal symbol of Lq for this choice of coordinates. Sufficient 
conditions for the uniqueness to the above Cauchy problem involve the properties of 
the roots r of the above polynomials r ^ Pq{-,T) as a function of ^. In the setting of 
redundant measurements, these conditions may hold for different values of q depending 
on the value of ^. This justifies the following definition. 

We assume the existence of a a finite covering {fij,} of the unit sphere (corre- 
sponding to 1^1 = 1) such that the following holds. For each z/, there exists q = q{v) 
and e > such that for each (x, t) close to and each ,^ G ^2;^, we have 

{i) T — )■ Pq{x, t, ^, t) has at most simple real roots and at most double complex roots 
(a) distinct roots ti and T2 satisfy |ri — > e > 
[in) non-real roots r satisfy |$5r| > e. 
Then we have the following result: 

Theorem A.l [Calderon's result for redundant systems.] Assume that N is non char- 
acteristic for the operators Lq at the origin and that for a finite covering {^lu} of the 
unit sphere, (i)-(ii)-(iii) above are satisfied. Then (148!) implies that u = in a full 
neighborhood ofO. 

Proof. The proof very closely follows that of Theorem 5 in [1D| . Note that the proof 
there requires all coefficients to be of class C°°. It has been shown in [25] (see also [SDl 
Notes of Chapter 2]) that such results also hold under less restrictions of smoothness. 
We do not track such results and assume that the coefficients are sufficiently smooth so 
that the technology developed in [10] applies. 

The differences between the proof of our theorem and that in [ID] are mostly nota- 
tional. To present these differences in a somewhat self-contained manner, we need to 
introduce most of the notation in [ID], to which we refer the reader for details that are 
not reproduced here. 

Let u{x,t) be of class with support contained in |x| < r and < t < T. Define 
w = e'^^^"*^ . Then there is a constant C independent of u such that for r, T, and k^^ 
sufficiently small, we have the Carleman-type inequality 

Q 

|a|<m 9=1 
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where we have defined 



\u\ 



T 

2 



2=/ w{t)\\u\\''dt, w{t) = e''^''-'^ , (50) 







with II • II the norm in the variable x. It is shown in [10|, p. 31] that the theorem 
follows from the above inequality, which we now prove. 

Since N is not characteristic for pg, the leading term of Pg in r"* does not vanish and 
we may divide Pg by such a term to consider leading symbols (still called Pg) of the form 

m 

with Qgj homogeneous polynomials in ^ of degree j. 

Let A = (1 — Aa;)2 be the operator with symbol (1 + |^P)^. For u with (small) 
compact support, we define 

Uj = A'^-Wt^, l<j<m. 
Let then U = (mi, . . . , UmY- We may recast LgU = fg as 

DtUj — Auj+i = 0, j < m 

m 

DtUm + '^Qgj{x,t,D^)A^~^Ura-j+l + Rgll = fg, 

i=i 

for Rg a pseudo-differential of order in x while Qqj{x, t, D^)A^~^ is a pseudo-differential 
of order 1 in x. The leading symbol of the above system for each q is therefore a m x m 
matrix rJ + hg with hg homogeneous of degree 1 in ^. We denote by Hq{t,x, D^) the 
pseudo-differential operator with symbol hg{t, x,^). Note that Det(rJ + /iq) = pg so that 
the roots of the system rJ + hg are those of Pg. 

From the derivation in |10l p. 33 & P-35], f H9|) follows (for and sufficiently 
small) from 

Q T 

<Cik-' + T')y2\\\iDt + Hg)U\\\^ + C [ wit)\\U\t,dt. (51) 

q=l Jo 

It thus remains to analyze the terms \\{Dt+Hq)U\\. The main idea of Calderon in [T7t [T8] 
is to diagonalize hq and analyze the corresponding diagonalized first-order system. For 
^ unit vectors close to some .^o, we find a non-singular smooth m x m matrix rq{x,t,C,) 
such that for all {x, t) close to 0, 

"^qhgfq Jg, 

with Jq in Jordan form. From the assumption (i), we obtain that for q = qiiy), then Jq 
is either a 1 x 1 matrix X{x,t,C,) or a 2 x 2 matrix with X{x,t,C,) on the diagonal, on 
the (2,1) entry, and 1 on the (1,2) entry. Once is defined for |,^| = 1, we extend it 
smoothly to all ^ as a homogeneous function of degree 0. This may be done locally on 
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cones in 1^1 > 1 but not necessarily globally. Following [101 p.34], the global obstructions 
may be overcome by replacing our finite covering {ft^} by a finer (sub) covering, which 
we still call {Qu} such that for each ^ G flu, there is a non-singular smooth matrix rg^ 
such that 

Tqyhqr"^ is in Jordan canonical form. 

We then define 0^(0 a C°° partition of unity subordinate to VL^ such that Yliy 4>l{0 = 
1 for 1^1 = 1 and such that each 0^ is extended to |^| > 1 smoothly as a homogeneous 
function of degree and such that = 1 1^1 > 1- Finally, we define as 

the corresponding pseudo-differential operator of order and := $^?7. 

Since [^u,Hq] is an operator of order 0, we verify as in [101 p. 35] that (1511) is a 
consequence of 

\\\Uu\\\^<C{k~'+T^)\\\{Dt + Hq^u))Uu\\\' + C r w{t)\\U\Wdt, (52) 

Jo 

which is to hold for each u. Note that the summation over q in (15T!) is a clear consequence 
of us obtaining the above inequality for q = q{v). For the rest of the derivation, q = q{v). 

Now, HqUy is unchanged if r^^ and jqy are changed outside of the support of 0^. 
We thus construct symbols that agree with hq on VL^ but that are globally reducible 
to a Jordan canonical form. Each is then extended to all ^ by a smooth function 
still called hy that is homogeneous of degree 1 for |^| > 1. We thus have globally 
defined symbols such that r^/i^r"^ = j^. We denote by H^, Ju, Ru, and the pseudo- 
differential operators of order 1, 1, 0, 0, respectively, corresponding to the symbols h^, 
ju, ru, and r'^ 

We finally define Vu = RuUu- The rest of the proof is then as in [10]. Following [101 
p. 35], we deduce that (152|) is a consequence of 

IIIKIII' < C{k-' + T^)\\m + (53) 

Now Ju obtained from pqi^u) for ^ G has diagonal entries of the form A{t) + iB{t) 
with A and B pseudo-differential operators of order 1 and B elliptic or vanishing when 
Ji, is a 1 X 1 matrix. Thus (!53|) is a consequence of [101 Lemma 2] and this concludes 
the proof of the theorem. □ 

A. 2 Redundant system for vector-valued functions 

The proof of the above theorem provides a local control encapsulated in fH9|) . Above, u 
is a scalar function. Such results may be extended to systems when u is vector-valued 
provided that the leading term of the system is diagonal. Indeed, ( 149|) may then be 
used component by component to provide a uniqueness result for the system. We now 
consider such an extension. 

Consider operators Lqr for 1 < q < Q and 1 < r < R and the system 

LqrUr = in fl > 0}, dl^+lUr = OB. V (1 {Xn+l = 0}, (54) 

with summation over repeated indexes and holding for all 1 < g < Q in the first equation 
and for all 1 < j < m — 1 and all 1 < r < i? in the second equation. Here V is still a 
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neighborhood of 0. We still call t = Xn+i and x = (xi, . . . , Xn) in an appropriate system 
with as before. We denote by Pgr = Pqr{x, t, ^, r) the principal symbol of Lgr for this 
choice of coordinates. As above, this uniqueness problem may be recast as: Assume 
that 

LqrUr = IB. V foT 1 < q < Q and = in {xn+i < 0}, (55) 

for 1 < r < R. Then = in a full neighborhood of 0. 

We sort the equations for 1 < g < Q such that the following holds. We assume that 
{1, . . . , Q} is decomposed into R connected components Ir of the form /(r — 1) + 1 < 
q < l{r) for a strictly increasing function /(r) from {1, . . . ,R} to N* with /(O) := — 1 
and 1{R) = Q. Then we assume that Pgr is of order m for q & Ir and of order at most 
m — 1 otherwise. When Q = R, then pg,. forms a determined systems with diagonal 
entries of order m and extra-diagonal entries of order at most m — 1. More generally, 
the matrix pgr is composed of operators of order m that are block- diagonal with blocks 
of size Ir X 1 and of other operators of order at most m — 1. 

Then, as we did before, we assume the existence of a finite covering {Qi,} of the unit 
sphere S"'"^ such that the following holds. For each u and r, there exists a g = g(z/, r) G Ir 
such that for each {x, t) close to and each ,^ G fi^, we have 

(0 Pqr{x,t,^,T) has at most simple real roots and at most double complex roots 
(ii) distinct roots Ti and T2 satisfy |ri — > e > 
{Hi) non-real roots r satisfy |53r| > e. 

Then we have the following result: 

Theorem A. 2 [Calderons result for Redundant systems.] Assume that N is non char- 
acteristic for the operators Lgr for q G Ir at the origin and that for a finite covering 
{flu} of the unit sphere, above are satisfied. Then fHS]) implies that u = 

in a full neighborhood of 0. 

Proof. The proof is similar to that of the preceding theorem. It is based on the 
following Carleman estimate 

R Q R 

^ ^ <C(A;-i+T2)^^|||Vn,||p, (56) 

''=1 |Q|<m q=l r=l 

whose proof is the same as that of (H^ . We leave the details to the reader. □ 

A. 3 Redundant systems in upper-triangular form 

In some applications, the leading term in the system may not be diagonal (or block- 
diagonal) but rather upper-triangular. In such a setting, unique continuation properties 
may still be valid provided that (a sufficiently large number of) the diagonal operators 
are elliptic. However, the corresponding complex roots may not longer be double. We 
have the following result, which as above, is stated for possibly redundant systems: 
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Theorem A. 3 Assume that N is non characteristic for the operators Lq at the origin 
and that there exists a finite covering {^y} of the unit sphere such that for each v, there 
exists q = q{v) such that 

Pq{x,t,^,T) has at most simple complex roots r satisfying |3r| > e 

for some e > 0. Then there exists a constant C independent of u such that for and 
T sufficiently small, we have 

Q 

\a\<m 9=1 

Proof. The main difference with respect to the preceding theorem is that we have a 
fuU control of the mth-order derivatives of u. We foUow the proof of Theorem lA. II until 
the very last step. The additional control of D'^u for |a| = m comes from the control 
of AUi, and DtU^, or equivalently, that of AVi, and DtV^. Now, hq(^y) is diagonalized 
with simple eigenvalues A(x, t,^) = — a(x, t,^) — ib{x,t,^). By assumption b{x,t,^) is 
uniformly bounded away from 0. We may then apply P0| (7.12)] to obtain that 



IIIAKlir + lllAKlir < C(l + kT')\\\{Dt + J.)K|ir. 



The proof of ( 1571) follows. □ 

As an application of the preceding result, let us consider the following corollary: 
Theorem A. 4 Consider the redundant system of equations 

Li Lq\ I Ui\ 

\ = mVn {xn+1 > 0}, a^+iMfe = onVn {xn+1 = 0}, (58) 

^Lg L2J \u2j 

for 1 < j <m — \ and 1 < k < 2, with the following assumptions. The operators Li and 
Lq are (vector-valued) Qi x 1 operators of order m, where Li satisfies the hypotheses of 
Theorem \A.l\ with Q = Qi. The operators L2 and L3 are Q2 x 1 operators of order m 
and at most m — 1, respectively. Moreover, L2 satisfies the hypotheses of Theorem \A.3\ 
with Q = Q2. Then [ui, U2) = in a full neighborhood of 0. 
The same result extends systems of the form 

(^ii)i<jj</jM = m V n > 0}, dl^^Uk = on Vn{xn+i = 0}, (59) 

for 1 < k < R, where Lij is (a vector-valued operator) of order m — 1 for i > j , La is 
(a vector-valued operator) of order m that satisfies the hypotheses of Theorem \A.1\ with 
an appropriate value of Q when all Lik are of order m — 1 for k ^ i and the hypotheses 
of Theorem \A . 31 with an appropriate value of Q when at least one operator Lik for k < i 
is of order m. 

Proof. First of all, since N is not characteristic for Li and L2, we observe that all 
derivatives of Uk vanish on {xn+i} so that both functions can be extended by for 
{xn+i < 0}. The above system of PDEs is thus satisfied on V with Uk = for Xn+i < 
and we wish to show that Uk = for Xn+i > 0. 
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The proof is based on a modification of the method showing Theorem I A. 1 1 from 
Indeed, we claim that for (mi,M2) with support contained in |x| < r and < t < T, 
then the theorem follows from 

2 

J2 $^|||/^"nfc|||2<C(A;"i + T2)|||Ln|||^ (60) 

a|<m k=l 

where we denote by Lu the left-hand side in f l58p . The proof is a straightforward 
modification of [101 P-31], which we reproduce here. Let ( he a. smooth function defined 
for t > equal to 1 for t < 2T/3 and for t > T. Let Uk = Qvk for v = {vi, V2) solution 
of Lv = 0. Then from f l60|) . using only the estimate with |a| = 0, we obtain that 



2T 



[ ' \\vfwdt<C{k~^ +T^) [ \\L{Cv)fwdt<C'{k~^ + T^)[ wdt 
Jo 

with C independent of k. This shows that 

T 



Jo 



and as /c — )■ cxD that J^^ WvW^dt — 0. 

It remains to show (l60l) . Using Theorem lA.ll we obtain that the left-hand-side is 
bounded by 

C{k-' + T^){\\\L^m\\\' + \\\L2U2\\\') <Cik~' + T^){\\\Lu\\\^ + \\\LoU2\\\' + \\\L,m\\\'). 
Now from Theorem IA.31 we obtain that 

\\\LoU2\\\^ <C{k~^ +T^)\\\L2U2\\\^. 

Moreover, we obtain that 

\\\Lsu^\\\' < C J2 \\\D''u^\\\'. 

\a\<m 

Therefore, for (/c~^ + T^) sufficiently small, we deduce ( l6Ql) . The generalization to (l59i) 
is mostly notational and is left to the reader. □ 
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